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Abstract 

The Effective One Body (EOB) formalism is an analytical approach 
which aims at providing an accurate description of the motion and radi- 
ation of coalescing binary black holes. We present a brief review of the 
basic elements of this approach. 



1 Introduction 

A network of ground-based interferometric gravitational wave (GW) detectors 
(LIGO/VIRGO/GEO/. . .) is currently taking data near its planned sensitivity. 
Coalescing black hole binaries are among the most promising, and most exciting, 
GW sources for these detectors. In order to successfully detect GWs from 
coalescing black hole binaries, and to be able to reliably measure the physical 
parameters of the source (masses, spins, . . .), it is necessary to know in advance 
the shape of the GW signals emitted by inspiralling and merging black holes. 
Indeed, the detection and subsequent data analysis of GW signals is made by 
using a large bank of templates that accurately represent the GW waveforms 
emitted by the source. 

Here, we shall introduce the reader to one promising strategy toward having 
an accurate analytical^ description of the motion and radiation of binary black 
holes, which covers all its stages (inspiral, plunge, merger and ring-down): the 
Effective One Body approach [J HI [31 g] . As early as 2000 |2] this method 
made several quantitative and qualitative predictions concerning the dynamics 
of the coalescence, and the corresponding GW radiation, notably: (i) a blurred 



^Here we use the adjective 'analytical' for methods that solve explicit (analytically given) 
ordinary differential equations (ODE), even if one uses standard (Rungo-Kutta-type) numeri- 
cal tools to solve them. The important point is that, contrary to 'id numerical relativity simula- 
tions, numerically solving ODE's is extremely fast, and can therefore be done (possibly even in 
real time) for a dense sampling of theoretical parameters, such as orbital (v = mi m^lM, . . .) 
or spin (a\ = Si/ Grn^ ,9i,ipi,...) parameters. 
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transition from inspiral to a 'plunge' that is just a smooth continuation of the 
inspiral, (ii) a sharp transition, around the merger of the black holes, between 
a continued inspiral and a ring-down signal, and (iii) estimates of the radiated 
energy and of the spin of the final black hole. In addition, the effects of the 
individual spins of the black holes were investigated within the EOB [S] and 
were shown to lead to a larger energy release for spins parallel to the orbital 
angular momentum, and to a dimensionless rotation parameter J / E"^ always 
smaller than unity at the end of the inspiral (so that a Kerr black hole can 
form right after the inspiral phase). All those predictions have been broadly 
confirmed by the results of the recent numerical simulations performed by several 
independent groups [ni[71[Hl[Sl[ini[ni[Il[Il[Tl[Tl[Tni[n] (for a review of 
numerical relativity results see [H]). Note that, in spite of the high computer 
power used in these simulations, the calculation of one sufficiently long waveform 
(corresponding to specific values of the many continuous parameters describing 
the two arbitrary masses, the initial spin vectors, and other initial data) takes on 
the order of two weeks. This is a very strong argument for developing analytical 
models of waveforms. 

Those recent breakthroughs in numerical relativity (NR) open the possibility 
of comparing in detail the EOB description to NR results. This EOB/NR com- 
parison has been recently initiated in several works [191 HOI HTJ [22l [23l [211 US . 
The level of analytical/numerical agreement is unprecedented, compared to 
what has been previously achieved when comparing other types of analytical 
waveforms to numerical ones. For instance, [24] found that the maximal de- 
phasing between a recent, very accurate NR simulation of 30 GW cycles during 
late inspiral [17] and EOB could be reduced below 10^"^ GW cycles. The same 
comparison exhibited also an excellent agreement between the amplitudes of 
the NR and EOB waveforms. If the reader wishes to put the EOB results in 
contrast with other (Post-Newtonian or hybrid) approaches he can consult, e.g.. 

Before reviewing some of the technical aspects of the EOB method, let us 
indicate some of the historical roots of this method. First, we note that the 
EOB approach comprises three, rather separate, ingredients: 

(1) a description of the conservative (Hamiltonian) part of the dynamics of 
two black holes; 

(2) an expression for the radiation-reaction part of the dynamics; 

(3) a description of the GW waveform emitted by a coalescing binary system. 

For each one of these ingredients, the essential inputs that are used in EOB 
developments are high-order post-Newtonian (PN) expanded results which have 
been obtained by many years of work, by many researchers (see references be- 
low). However, one of the key ideas in the EOB philosophy is to avoid using 
PN results in their original 'Taylor-expanded' form (i.e. cq + ci u + C2 f ^ -I- 
C3 z)^ -|- • • • -|- c„ w"), but to use them instead in some resummed form (i.e. some 
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non-polynomial function of w, defined so as to incorporate some of the expected 
non-perturbative features of the exact result). The basic ideas and techniques 
for resumming each ingredient of the EOB are different and have different his- 
torical roots. Concerning the ingredient (1), i.e. the EOB Hamiltonian, it was 
inspired by an approach to electromagnetically interacting quantum two-body 
systems introduced by Brezin, Itzykson and Zinn- Justin [21] . The resummation 
of the ingredient (2), i.e. the EOB radiation- reaction force was inspired by 
the Pade resummation of the flux function introduced by Damour, Iyer and 
Sathyaprakash [52] ■ As for the ingredient (3), i.e. the EOB description of the 
waveform emitted by a coalescing black hole binary, it was mainly inspired by 
the work of Davis, RufSni and Tiomno |35] which discovered the transition be- 
tween the plunge signal and a ringing tail when a particle falls into a black hole. 
Additional motivation for the EOB treatment of the transition from plunge to 
ring-down came from work on the, so-called, 'close limit approximation' j34|. 

Let us finally note that the EOB approach has been recently improved by 
following a methodology consisting of studying, element by element, the physics 
behind each feature of the waveform, and on systematically comparing vari- 
ous EOB-based waveforms with 'exact' waveforms obtained by NR approaches. 
Among these 'exact' NR waveforms, it has been useful to consider the small- 
mass-ratio limit ly = mi m2/(mi +1112)'^ ^ 1, in which one can use the well con- 
trollable 'laboratory' of numerical simulations of test particles (with an added 
radiation-reaction force) moving in black hole backgrounds [23) . 

2 Motion and radiation of binary black holes: 
post-Newtonian expanded results 

Before discussing the various resummation techniques used in the EOB ap- 
proach, let us briefly recall the 'Taylor-expanded' results that have been ob- 
tained by pushing to high accuracies the post-Newtonian (PN) methods. 

Concerning the orbital dynamics of compact binaries, we recall that the 
2.5PN-accurat41 equations of motion have been derived in the 1980's [SS] [Ml 
[37l|38]. Pushing the accuracy of the equations of motion to the 3PN (~ (v/c)^) 
level proved to be a non-trivial task. At first, the representation of black holes 
by delta-function sources and the use of the (non diffeomorphism invariant) 
Hadamard regularization method led to ambiguities in the computation of the 
badly divergent integrals that enter the 3PN equations of motion [3S1 [3D] . This 
problem was solved by using the (diffeomorphism invariant) dimensional regu- 
larization method (i.e. analytic continuation in the dimension of space d) which 
allowed one to complete the determination of the 3PN-level equations of mo- 
tion [m [42] . They have also been derived by an Einstein-Infeld-Hoffmann-type 
surface- integral approach [43] . The 3.5PN terms in the equations of motion are 
also known [31[^[3S]. 

^As usual 'n-PN accuracy' means that a result has been derived up to (and including) 
terms which are ~ (v/c)-^" ~ (GM/c^r)^ fractionally smaller than the leading contribution. 
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Concerning the emission of gravitational radiation, two different gravitational- 
wave generation formalisms have been developed up to a high PN accuracy: 
(i) the Blanchet-Damour-Iyer formalism gTl IIHI SHI [SUl [5T1 [SH combines 
a multipolar post-Minkowskian (MPM) expansion in the exterior zone with a 
post-Newtonian expansion in the near zone; while (ii) the Will-Wiseman-Pati 
formalism [SH [511 ESI Si] uses a direct integration of the relaxed Einstein equa- 
tions. These formalisms were used to compute increasingly accurate estimates of 
the gravitational waveforms emitted by inspiralling binaries. These estimates in- 
clude both normal, near-zone generated post-Newtonian effects (at the IPN [15] . 
2PN [571 [551 [M]) and 3PN [5S1[SD| levels), and more subtle, wave-zone generated 
(Hnear and non-Hnear) 'tail effects' [5T1 [5T1 [SH [53] . However, technical prob- 
lems arose at the 3PN level. The representation of black holes by 'delta- function' 
sources causes the appearance of dangerously divergent integrals in the 3PN mul- 
tipole moments. The use of Hadamard (partie finie) regularization did not allow 
one to unambiguously compute the needed 3PN-accurate quadrupole moment. 
Only the use of the (formally) diffeomorphism-invariant dimensional regular- 
ization method allowed one to complete the 3PN-level gravitational-radiation 
formalism [63] . 

The works mentioned in this Section (see for a detailed account and more 
references) finally lead to PN-expanded results for the motion and radiation of 
binary black holes. For instance, the equations of motion are given in the form 
(a= 1,2; z = 1,2,3) 

^ X icons . A'lRB. f-tS 

where 

A™"*^ = Ao + c-^A2 + c-^Ai + c-^Ae , (2) 
denotes the 'conservative' 3PN-accurate terms, while 

denotes the time-asymmetric contibutions, linked to 'radiation reaction'. 

On the other hand, if we consider for simplicity the inspiralling motion of 
a quasi-circular binary system, the essential quantity describing the emitted 
gravitational waveform is the phase (p of the quadrupolar gravitational wave 
amplitude h{t) ~ a{t) cos{(j){t) -\- 6). PN theory allows one to derive several 
different functional expressions for the gravitational wave phase 0, as a function 
either of time or of the instantaneous frequency. For instance, as a function 
of time, (j) admits the following explicit expansion in powers of = vc^ltc — 
t)/5GM (where tc denotes a formal 'time of coalescence', M = mi -\- m2 and 
u = mi m2/M^) 

cj,{t) = 0, - v-^ 05/8 + Y^^a^ + In e) j , (4) 

with some numerical coefficients a„, aj^ which depend only on the dimensionless 
(symmetric) mass ratio v = mi m2/M'^ . The derivation of the 3.5PN-accurate 
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expansion ^ uses both the 3PN-accurate conservative acceleration ^ and a 
3.5PN extension of the (fractionally) IPN-accurate radiation reaction acceler- 
ation ([3]) obtained by assuming a balance between the energy of the binary 
system and the gravitational- wave energy flux at infinity (see, e.g., [64] ) . 

3 Conservative dynamics of binary black holes: 
the Effective One Body approach 

The PN-expanded results briefly reviewed in the previous Section are expected 
to yield accurate descriptions of the motion and radiation of binary black holes 
only during their early inspiralling stage, i.e. as long as the PN expansion pa- 
rameter 7e — GM/c^R (where R is the distance between the two black holes) 
stays signiflcantly smaller than the value ^ i where the orbital motion is ex- 
pected to become dynamically unstable ('last stable circular orbit' and begin- 
ning of a 'plunge' leading to the merger of the two black holes). One needs a 
better description of the motion and radiation to describe the late inspiral (say 
7e > J2), as well as the subsequent plunge and merger. One possible strat- 
egy for having a complete description of the motion and radiation of binary 
black holes, covering all the stages (inspiral, plunge, merger, ring-down), would 
then be to try to 'stitch together' PN-expanded analytical results describing the 
early inspiral phase with 3d numerical results describing the end of the inspi- 
ral, the plunge, the merger and the ring-down of the flnal black hole, see, e.g., 
Refs. [651 [20]. 

However, we wish to argue that the EOB approach makes a better use of all 
the analytical information contained in the PN-expanded results ([l])-©. The 
basic claim (first made in [Il[2) is that the use of suitable resummation methods 
should allow one to describe, by analytical tools, a sufficiently accurate approx- 
imation of the entire waveform, from inspiral to ring-down, including the non- 
perturbative plunge and merger phases. To reach such a goal, one needs to make 
use of several tools: (i) resummation methods, (ii) exploitation of the flexibility 
of analytical approaches, (iii) extraction of the non-perturbative information 
contained in various numerical simulations, (iv) qualitative understanding of 
the basic physical features which determine the waveform. 

Let us start by discussing the first tool used in the EOB approach: the 
systematic use of resummation methods. Two such methods have been employed 
(and combined), and some evidence has been given that they do significantly 
improve the convergence properties of PN expansions. The first method is the 
use of Fade approximants. It has been shown in Ref. [32j that near-diagonal 
Pade approximants of the radiation reaction forc^ T seemed to provide a good 



^We henceforth denote by the Hamiltonian version of the radiation reaction term A^^, 
Eq. Jsj, in the (PN-expanded) equations of motion. It can be heuristically computed up to 
(absolute) 5.5PN |59l 1661 and even 6PN |67l order by assuming that the energy radiated 
in gravitational waves at infinity is balanced by a loss of the dynamical energy of the binary 
system. 
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representation of T down to the last stable orbit (which is expected to occur 
when R ~ ^GMjc? , i.e. when 7e ~ i). The second method is a novel approach 
to the dynamics of compact binaries, which constitutes the core of the Effective 
One Body (EOB) method. 

For simplicity of exposition, let us first explain the EOB method at the 
2PN level. The starting point of the method is the 2PN-accurate Hamilto- 
nian describing (in Arnowitt-Deser-Misner-type coordinates) the conservative, 
or time symmetric, part of the equations of motion ([1]) (i.e. the truncation 
= ^0 + 0-2^2 + c-4^4 of Eq. dm) say i/2PN('Zi - 92,^1,^2). By going 
to the center of mass of the system (pi + P2 = 0), one obtains a PN-expanded 
Hamiltonian describing the relative motion, q = Qi q2, P = Pi = —P2- 

Hgt'"'iQ,P) = HoiQ,P) + \H2{q,p) + l^Hi{q,p) , (5) 

where Ho(q,p) = J- p'^ + ^MJi. (with M = mi + m2 and ^ = mim2/M) 
corresponds to the Newtonian approximation to the relative motion, while H2 
describes IPN corrections and H4 2PN ones. It is well known that, at the New- 
tonian approximation, HQ(q,p) can be thought of as describing a 'test particle' 
of mass fi orbiting around an 'external mass' GM. The EOB approach is a gen- 
eral relativistic generalization of this fact. It consists in looking for an 'external 
spacetime geometry' g^ (x"^ ; GM) such that the geodesic dynamics of a 'test 
particle' of mass /i within g^{x^, GM) is equivalent (when expanded in powers 
of 1/c^) to the original, relative PN-expanded dynamics ([5]). 

Let us explain the idea, proposed in [1], for establishing a 'dictionary' be- 
tween the real relative- motion dynamics, ([5]), and the dynamics of an 'effective' 
particle of mass n moving in g'^{x'^, GM). The idea consists in 'thinking quan- 
tum mechanically 0. Instead of thinking in terms of a classical Hamiltonian, 
H{q,p) (such as TJjpn'™' ^Iq. ([5])), and of its classical bound orbits, we can 
think in terms of the quantized energy levels E{n, £) of the quantum bound 
states of the Hamiltonian operator H{q,p). These energy levels will depend on 
two (integer valued) quantum numbers n and £. Here (for a spherically sym- 
metric interaction, as appropriate to //relative -j^ £ parametrizes the total orbital 
angular momentum (L^ = £{£ + 1) ti^), while n represents the 'principal quan- 
tum number' n = £ + rir + 1, where Ur (the 'radial quantum number') denotes 
the number of nodes in the radial wave function. The third 'magnetic quantum 
number' m (with —£< m < £) does not enter the energy levels because of the 
spherical symmetry of the two-body interaction (in the center of of mass frame). 
For instance, a non-relativistic Coulomb (or Newton!) interaction 

^o = ^P^ + ^ (6) 
2^ \q\ 



*This is related to an idea emphasized many times by John Archibald Wheeler: quantum 
mechanics can often help us in going to the essence of classical mechanics. 
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gives rise to the well-known result 



which depends only on n (this is the famous Coulomb degeneracy). When 
considering the PN corrections to Hq, as in Eq. ([5]), one gets a more complicated 
expression of the form 
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^2PN ("'^'•-^2^^ 



1 + (211 + '^20 \ 

\ n£ Ti? / 

/ Cl3_ C22 C3]_ c^N 



(8) 



where we have set a = GMfj,/h = Gmim2/fi, and where we consider, for 
simplicity, the (quasi-classical) limit where n and £ are large numbers. The 
2PN-accurate result ([S]) had been derived by Damour and Schafer ^S] as early 
as 1988. The dimensionless coefficients Cpq are functions of the symmetric mass 
ratio v = n/M , for instance C40 = |(145— I'bv + v^). In classical mechanics [i.e. 
for large n and 1), it is called the 'Delaunay Hamiltonian', i.e. the Hamiltonian 
expressed in terms of the action variable)^ J ^ £fi = ^ § p^p dip, and N = nh^ 
Ir + J, with Ir = § Pr dr. 

The energy levels ([5]) encode, in a gauge-invariant way, the 2PN-accurate 
relative dynamics of a 'real' binary. Let us now consider an auxiliary problem: 
the 'effective' dynamics of one body, of mass fi, following a geodesic in some 
'external' (spherically symmetric) metric!^ 

dx^' dx" -A{R) dT^ + B{R) dR^ + R^{de^ + sin^ 6 dip^) . (9) 

Here, the a priori unknown metric functions A{R) and B{R) will be constructed 
in the form of an expansion in GM/c^R: 

GM fGM\^ /GmV 



GM (GM\ 
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BiR) = l + ^i^+^2^^j (10) 

where the dimensionless coefficients a„,6„ depend on v. From the Newtonian 
limit, it is clear that we should set ai — —2. By solving (by separation of 
variables) the 'effective' Hamilton-Jacobi equation 



dScS dScS ,22 n 

=°' 



^We consider, for simplicity, 'equatorial' motions with m = £, i.e., classically, 9 = ^. 

®It is convenient to write the 'external metric' in Schwarzschild-like coordinates. Note that 
the external radial coordinate _R differs from the two-body ADM-coordinate relative distance 
^ADM = |q|. The transformation between the two coordinate systems has been determined 
in Refs. 
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ScS ^ -ScS t + Jeff ip + Scs{R) , (11) 

one can straightforwardly compute (in the quasi-classical, large quantum num- 
bers limit) the Delaunay Hamiltonian £off(-^cff, Jeff), with N^ff = ricff K Jeff = 
4fffi (where A^eff - Jeff + if, with if = ^ § pf dR, Pf = dSeffiR)/dR). 
This yields a result of the form 



fofT(ricff , 4fr) = /^c^ - - M 



2 «e'ff 



„cff 


„ofI \ 
"off/ 


'T'cff^ff 


„ofT 


„cfT 
1 ''22 


nofT^eff 


7)2 «2 



(12) 



-I- 



„cff „off 
•^Sl I '-40 



"efT^oft 



where the dimensionless coefficients c^^ are now functions of the unknown co- 



efficients a„,6„ entering the looked for 'external' metric coefficients (|TO)l . 

At this stage, one needs (as in the famous AdS/CFT correspondence) to 
define a 'dictionary' between the real (relative) two-body dynamics, summarized 
in Eq. ([8]), and the effective one-body one, summarized in Eq. p2|) . As, on both 
sides, quantum mechanics tells us that the action variables are quantized in 
integers (A^reai = nh, Neff = Ueffh, etc.) it is most natural to identify n = Ueff 
and £ = ieff. One then still needs a rule for relating the two different energies 
-^rcaj*'™ ^cff- n proposed to look for a general map between the 

real energy levels and the effective ones (which, as seen when comparing ([8]) 
and (|12p . cannot be directly identified because they do not include the same 
rest- mass contributiorQ) , namely 



C „ / TTTclativc \ TT^rclativc / TT^rclativc / TT^rclativc ^ ^ 

= / ( Ei^] ^ 1 + + a, ( 

(13) 

The 'correspondence' between the real and effective energy levels is illustrated 
in Fig. 1. 

Finally, identifying £efi{n,t)/^(? to f [Elf^^^° / ^xc^) yields six equations, re- 
lating the six coefficients Cp^(a2, 03; 6i, 62) to the six Cpq{v) and to the two 
energy coefficients ai and a2- It is natural to set 61 = +2 (so that the lin- 
earized effective metric coincides with the linearized Schwarzschild metric with 
mass M = mi -f 7712). One then finds that there exists a unique solution for the 
remaining five unknown coefficients 02,03,62,^1 and a2- This solution is very 
simple: 

02 = , 03 = 2i/ , 62 = 4 — 61/ , ai = — , a2 = . (14) 
Note, in particular, that the map between the two energies is simply 

c TT^iclativc / TT^rclativc \ „ ,m2 „4 ™2 „4 



■^Indeed Ej°*f' = A/c^ + B^iJ^f = Mc^ -|- Newtonian terms + IPN/c^ -I- • ■ ■ , while 
^^offoctivc = + AT + lPN/c2 + ■ . ■ . 
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£ = f{E) 




n + 1,^ + 1 



Figure 1: Sketch of the correspondence between the quantized energy levels of 
the real and effective conservative dynamics, n denotes the 'principal quantum 
number' (n = + £+ 1, with Ur = 0, 1, . . . denoting the number of nodes in the 
radial function), while £ denotes the (relative) orbital angular momentum {L^ = 
£{t+l) It?). Though the EOB method is purely classical, it is conceptually useful 
to think in terms of the underlying (Bohr-Sommerfeld) quantization conditions 
of the action variables and J to motivate the identification between n and (. 
in the two dynamics. 

where s = (^roli)^ = (^^^^ + ^^Jeif is Mandelstam's invariant. Note also 
that, at 2PN accuracy, the crucial '^go*' metric coefficient A{R) (which fully 
encodes the energetics of circular orbits) is given by the remarkably simple PN 
expansion 

^2PN(i?) = l-2M + 2^y^i^ (16) 

where u = GM/c^R and v = f^/M = nii 7712/ (mi + TO2)^- 

The dimensionless parameter v = fi/M varies between (in the test mass 
limit mi ^ 111,2) and j (in the equal-mass case mi = m2). When 1/ ^ 0, 
Eq. (jl6[) yields back, as expected, the well-known Schwarzschild time-time met- 
ric coefficient -ffoo''" = 1 - 2u = 1 - 2GM /c^R. One therefore sees in Eq. (fT6|) 
the role of as a deformation parameter connecting a well-known test-mass 
result to a non trivial and new 2PN result. It is also to be noted that the 
IPN EOB result AiPN(i?) — l~2u happens to be i^-independent, and therefore 
identical to A^'^'^^ = 1 — 2u. This is remarkable in view of the many non-trivial 
i^-dependent terms in the IPN relative dynamics. The physically real IPN u- 
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dependence happens to be fully encoded in the function f{E) mapping the two 
energy spectra given in Eq. above. 

Let us emphasize the remarkable simplicity of the 2PN result (fT6)) . The 
2PN Hamiltonian ([5]) contains eleven rather complicated z/-dependent terms. 
After transformation to the EOB format, the dynamical information contained 
in these eleven coefficients gets compactified into the very simple additional 
contribution +2h'u^ in A{R), together with an equally simple contribution in 
the radial metric coefficient: { A{R) B{R)) 2Pti = 1 — 6^m^. This compactifi- 
cation process is even more drastic when one goes to the next (conservative) 
post-Newtonian order: the 3PN level, i.e. additional terms of order 0{l/c^) in 
the Hamiltonian ([5]). As mentioned above, the complete obtention of the 3PN 
dynamics has represented quite a theoretical challenge and the final, resulting 
Hamiltonian is quite complicated. Even after going to the center of mass frame, 
the 3PN additional contribution HQ{q,p) to Eq. ([5|) introduces eleven new 
complicated z^-dependent coefficients. After transformation to the EOB format 
[3], these eleven new coefficients get 'compactified' into only three additional 
terms: (i) an additional contribution to A(R), (ii) an additional contribution 
to B{R), and (iii) a 0{p'^) modification of the 'external' geodesic Hamiltonian. 
For instance, the crucial 3PN ^qq' metric coefficient becomes 

A^pj^iR) ^l-2u + 2uu^ + a4vu^ , (17) 

where 

94 41 

04 ^ - ^ TT^ ~ 18.6879027 . (18) 

Remarkably, it is found that the very simple 2PN energy map Eq. (|15p does not 
need to be modified at the 3PN level. 

The fact that the 3PN coefficient 04 in the crucial 'effective radial potential' 
^3Pn(-R), Eq. ^T7\ . is rather large and positive indicates that the i^-dependent 
nonlinear gravitational effects lead, for comparable masses {1/ i), to a last 
stable (circular) orbit (LSO) which has a higher frequency and a larger binding 
energy than what a naive scaling from the test-particle limit (ly — > 0) would 
suggest. Actually, the PN-expanded form (fT7|) of A3p-^{R) does not seem to 
be a good representation of the (unknown) exact function Aeob(J^) when the 
(Schwarzschild-like) relative coordinate R becomes smaller than about QGM/(? 
(which is the radius of the LSO in the test-mass limit). It was therefore sug- 
gested |3| to further resunjl yl3PN(i?) by replacing it by a suitable Pade (P) 
approximant. For instance, the replacement of A3PN(i?) by 

Aim-pi\A,Mm- ,^,J^:^^,^^, (19) 

ensures that the v ^ j case is smoothly connected with the v = {] limit. 



*The PN-expanded EOB building blocks A{R), B{R), . . . already represent a resummation 
of the PN dynamics in the sense that they have compactified the many terms of the original 
PN-expanded Hamiltonian within a very concise format. But one should not refrain to further 
resum the EOB building blocks themselves, if this is physically motivated. 
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The use of ([T^ was suggested before one had any (rehable) non-perturbative 
information on the binding of cfose black hole binaries. Later, a comparison 
with some 'waveless' numerical simulations of circular black hole binaries [69j 
has given some evidence that is physically adequate. In Refs. [H HH] it was 
also emphasized that, in principle, the comparison between numerical data and 
EOB-based predictions should allow one to determine the effect of the unknown 
higher PN contributions to Eq. (fT7| . For instance, one can add a 4PN-like 
term +051^^^ in Eq. (|17p . and then Fade the resulting radial function, say 
A\ = P|[^3PN + a^vu^]. Comparing the predictions of ^4(05] to numerical 
data might then determine what is the physically preferred 'effective' value of 
the unknown coefficient 05. This is an example of the useful flexibility [70j 
of analytical approaches: the fact that one can tap numerically-based, non- 
perturbative information to improve the EOB approach. 

4 Description of radiation-reaction effects in the 
Effective One Body approach 

In the previous Section we have described how the EOB method encodes the 
conservative part of the relative orbital dynamics into the dynamics of an 'ef- 
fective' particle. Let us now briefly discuss how to complete the EOB dynamics 
by defining a resummed expression describing radiation reaction effects. One 
is interested in circularized binaries, which have lost their initial eccentricity 
under the influence of radiation reaction. For such systems, it is enough to 
include a radiation reaction force in the equation of motion. More precisely, 
we are using phase space variables r,pr,(p,p^ associated to polar coordinates 
(in the equatorial plane 9 = Actually it is convenient to replace the radial 
momentum pr by the momentum conjugate to the 'tortoise' radial coordinate 
i?* = J dR{B/A)'^/^, i.e. Pr, = {A/Bf/^Pn. The real EOB Hamiltonian is 
obtained by first solving Eq. p5|l to get = \fs in terms of feff, and then 

by solving the effective Hamiltonian-Jacobi equatior|^ to get EeS in terms of the 
effective phase space coordinates qcff and Pcff- The result is given by two nested 
square roots (we henceforth set c— 1): 



with 2:3 = 2j/ (4 — 2>v) . Here, we are using suitably rescaled dimensionless (effec- 
tive) variables: r = R/GM, p^, = PR,/fJ; P<p = Pip/fJ- GM, as well as a rescaled 



Completed by the ©(p^) terms that must be introduced at 3PN. 




(20) 



where 




(21) 
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time t = T/GM. This leads to equations of motion of the form 



dr 
'di 

dpr, 

dt 



dH 



EOB 



dPr, 



1/2 



EOB 



Or 



^ _ dH^OB 
dt dp^ 



dp^ 

~dr 



(22) 



As explained above the EOB metric function A(r) is defined by Pade resumming 
the Taylor-expanded result pop obtained from the matching between the real 
and effective energy levels (one uses a similar Pade resumming for B{r), or 
rather for D{r) = A(r)B{r)). One similarly needs to resum the ip component 
of the radiation reaction which has been introduced on the r.h.s. of the last 
equation (|22p . During the quasi-circular inspiral J^^ is known (from the PN 
work mentioned in Section 2 above) in the form of a Taylor expansion of the 
form 

32 ^4 ^Taylor (23) 



-^Taylor 



where = VLr^ 



ril;^/^ (with -0 defined as in Eq. (22) of [TT]), and 



^Taylor = I + A2{v) + A^{v) + Ai{v) v'^ + A^{v) v"^ 

+ AQ{v,\ogv)v^ + A^{v)v^ + As{v = O,log^;)^;^ (24) 

where we have added to the known 3.5PN-accurate comparable-mass result the 
small-mass-ratio 4PN contribution [72] . 

Following [32] ) one resums ^""'"'^^1°'' by using the following Pade resummation 
approach. First, one chooses a certain number fpoio which is intended to repre- 
sent the value of at which the exact angular momentum flux would become 
infinite if one were to formally analytically continue along unstable circular 
orbits below the Last Stable Orbit (LSO): then, given Wpoic, one defines the 
resummed F{vtp) as 



F 



resummed 



'^^polc 



Vpolc 



(25) 



where denotes a (4,4) Pade approximant. It has recently been shown [M] 
that the flexibility in the choice of Wpoio could be advantageously used to get 
a close agreement with NR data. The quality of the agreement with NR data 
obtained by such a 'flexed' Pade resummation is illustrated (in the small mass- 
ratio limit) in Fig. 2. 
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Figure 2: This figure illustrates how the Pade resummation of the (Newton- 
normalized) 3PN-expanded angular momentum flux F"^^^^°'^{v) — F{v)/Fn{v) 
nicely agrees with the exact, numerical flux when using a suitably 'flexed' value 
of Vpo\e- The solid curve is the exact, numerical flux [73]. The empty circles 
denote the (3, 3) Fade resummation using Wpoic — 0.53. The dashed curve, below 
the other two, denote the (3,3) Fade resummation using the value Wpoio = l/'s/Sj 
as originally suggested [32] . 

Let us note that the EOB method is not restricted to planar motions. In 
particular, the entire method has been extended to the case of circularized 
motions of spinning black holes [H [SJ [7S]. In this case, one must work with 
more phase space variables, q, p, Si and S2, and include spin effects in the 
radiation reaction. 

5 Effective One Body waveforms 

To end this brief review, let us sketch the definition of accurate waveforms 
within the EOB approach. The construction of EOB waveforms is based on 
several separate inputs: 

• high-accuracy, FN-expanded inspiralling waveforms as basic inputs (see 
references above); 

• a specific resummation of the inspiral-plus-plunge waveform h™'^^^^'^^'^ (t) 
(which includes the resummation of an infinite number of 'leading log- 
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arithms' entering the Multipolar-post-Minkowskian expansion of tail ef- 
fects, the resummation of the 5PN-accurat43 'non-tail' i = 2, m — 2 
waveforrrF^. and the parametrization of some non-quasi-circular effects) 

(1^) = -8 /I Hr. F,, fg^^ ; (26) 

• a simplified representation of the transition between plunge and ring-down 
by smoothly matching, on a (2p -|- l)-toothed "comb" (t„i — pS, . . . ,t,n ~ 
5, tm,tm+S, . . . , tm+pS) Centered around a matching time tm, the inspiral- 
plus-plunge waveform to a ring-down waveform, made of the superposition 
of severao quasi-normal-mode complex frequencies, 

(^) = E e-«(*-*-) , (27) 

V / AT 

with (Tj^ = ajv + * ojn, and where the label N refers to indices {£, £' , to, n), 
with {i,m) = (2,2) being the Schwarzschild-background multipolarity of 
the considered (metric) waveform him, with n = 0, 1, 2 . . . being the 'over- 
tone number' of the considered Kerr-background Quasi-Normal-Mode, and 
£' the degree of its associated spheroidal harmonics Se'm{ao;9). As dis- 
cussed in |2] and [53] the physics of the transition between plunge and ring- 
down (which was first understood in the classic work of Davis, Ruffini and 
Tiomno [33j) suggests to choose as matching time tm, in the comparable- 
mass case, the EOB time when the EOB orbital frequency ^l{t) reaches 
its maximum value. 

Here, we have been describing the 3+^-PN-accurate resummed EOB wave- 
form introduced in [531 IMl ■ A less accurate, restricted (0-PN) EOB waveform 
(together with a simplified matching procedure; 6^0) has been used by Buo- 
nanno and collaborators [21]. 

Finally, one defines a complete, quasi-analytical EOB waveform (covering 
the full process from inspiral to ring-down) as: 

hf^^it) = eitm - t) 4"2^P'^"SC(^) ^ ^(^ _ ;,™gdow„(^^ ^ (28) 

where 9{t) denotes Heaviside's step function. 

An example of this complete EOB waveform is represented in Fig. 3. 



^^The 3PN-accurate part of the 'non-tail' waveform is known in the comparable-mass case 
V ^0. The 4PN and 5PN pieces are only known in the test- mass limit i/ — ► [721 173) . 

^^See I74| for an independent derivation of the non-resummed 3PN-accurate I = 2, m = 2 
waveform. 

^^Refs. 1231 1251 use p = 2, i.e. a 5-teethed comb, and, correspondingly, 5 positive-frequency 
Kerr Quasi-Normal Modes. 
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Figure 3: This figure illustrates a complete resummed EOB quadrupolar {i = 2, 
m = 2) waveform (|28|) (using = 20, aim = 30 and a suitably flexed Upoic), 
with, about, 29 GW cycles of inspiral, ~ 1 GW cycle during plunge, and ~ 4 
GW cycles of ring-down. Ref. [H] has shown that the inspiral part of the EOB 
waveform agreed with the Caltech-Cornell NR data [17] within ±0.001 GW 
cycles, while Ref. [25] has shown that the late-inspiral-plunge-ringdown part 
of the EOB waveform agrees with an Albert Einstein Institute NR waveform 
within ±0.005 GW cycles. 

Recent comparisons between EOB waveforms and full, 3-dimensional NR 
waveforms have exhibited an excellent level of agreement. After a prelimi- 
nary comparison done in [19j . Buonanno et al. [21j compared restricted EOB 
waveforms to NR waveforms computed by the NASA-Goddard group. In the 
equal- mass case — 1/4), they found that the dephasing between (restricted) 
EOB and NR waveforms (covering late inspiral, merger and ring-down) stayed 
within ±0.030 GW cycles over 14 GW cycles. In the case of a mass ratio 4 : 1 
(i/ = 0.16), the dephasing stayed within ±0.035 GW cycles over 9 GW cycles. 
The resummed EOB waveform has been compared to two different (equal-mass) 
NR waveforms: (i) in the comparison with the very accurate inspiralling simula- 
tion of the Caltech-Cornell group [T7] the dephasing stayed smaller than ±0.001 
GW cycles over 30 GW cycles (and the amplitudes agreed at the ~ 10"'^ level) 
|24j ; (ii) in the comparison with a late-inspiral-merger-ringdown NR waveform 
computed by the Albert Einstein Institute group, the dephasing stayed smaller 
than ±0.005 GW cycles over 12 GW cycles. See Ref. for another EOB/NR 
comparison. All these comparisons have used the natural flexibility of the EOB 
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formalism [H [7D] , i. e. the possibility to tune EOB parameters representing yet 
uncalculated effects (such as 05 in A(r) or Wpoie)- 

We hope to have shown here the usefulness of the EOB formalism in ac- 
curately describing the general relativistic motion and radiation of coalescing 
black hole binaries. 
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